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Abstract 

Properties of the domain wall (kink) solution in the 5 dimensional Randall- 
Sundrum model are examined both analytically and numerically. The con- 
figuration is derived by the bulk Higgs mechanism. We focus on 1) the 
convergence property of the solution, 2) the stableness of the solution, 3) the 
non-singular property of the Riemann curvature, 4) the behaviours of the 
warp factor and the Higgs field. It is found that the bulk curvature changes 
the sign around the surface of the wall. We also present some exact solu- 
tions for two simple cases: a) the no potential case, b) the cosmological term 
dominated case. Both solutions have the (naked) curvature singularity. We 
can regard the domain wall solution as a singularity resolution of the exact 
solutions. 
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1 Introduction 

There exist two standpoints in the treatment of the domain world physics. 
One is that the geometry should be singular and the domain is regarded as a 
defect. In the original work by Randall and Sundrum|l|, the walls stand on 
the fixed points of the S1/Z2 orbifold in the form of (^-function. Recently the 
renormalization of a bulk-boundary system is discussed in this standpoint [0 
where the conical singularity at the center of the extra 2 dim space is regarded 
as a domain of the vortex type. The other standpoint is that the geometry 
should be non- singular. The domain is regarded as a soliton. This approach 
looks natural if the domain world is "derived" from the more fundamental 
theory of the D-brane. Similar situations occurred in the past literature, such 
as the relation between the Dirac string || and 'tHooft-Polyakov string ||, ||]. 
(For the situation about the (topological) defect and the soliton in relation 
to the domain world, see a good review ||.) At present both standpoints 
look important to understand the brane world physics. We take here the 
latter standpoint. 

In ref.jTJ, the domain wall configuration of the RS-model is realized as 
a soliton (kink) solution in the bulk (5D) Higgs potential. It has some ad- 
vantageous points, compared with the (^-function description, such as non- 
singularity and stability. The solution is obtained in the form of the infinite 
power-series of some hyperbolic function. The convergence of the coefficient- 
series is crucial for the boundary condition (b.c.) to be satisfied. It was 
checked by explicitly calculating the coefficients at the 2nd order. We present 
here the 6th order calculation result, and reconfirm the convergence property 
further strongly. The main purpose of this paper is to strengthen the content 
of ref . J7| by presenting the various results in a concrete way. Some physical 
quantities, such as the bulk scalar curvature, are obtained. Very interest- 
ingly, the curvature changes its sign near the "surface" of the domain wall. 
In order to clarify the structure of the solution, we first present some exact 
solutions for simple cases. They clearly show the origin of some integration 
constants and free parameters. These exact solutions have (naked) singu- 
larities. They tell us the Higgs potential is important for the non-singular 
property of the configuration. It plays the role of singular resolution. 

We consider one- wall model which was considered in ||. An interesting 
stable (kink) solution exists for a family of vacua. As explained in [7], the 
solution does not miss the key points of the original one. Similar analysis 
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was successfully done in the 6 dim mo del [ 

As some recent related works, we find 

In Sec. 2, the RS domain wall model is explained. The exact solution 
is presented for the no-potential case in Sec. 3, and for the cosmological- 
constant dominated case in Sec. 4. The general case of the Higgs potential 
is examined in Sec. 5, where the domain- wall solution is obtained as a one- 
parameter family of kink solutions. It is the analytical solution. The concrete 
values of parameters and coefficients are obtained in the 6th order calculation. 
Numerical analysis also confirms the obtained solution. 

2 Randall- Sundrum model with the bulk Higgs 
field 

We consider, as the brane world, the following 5D gravity-Higgs theory. 

3[GabM = J d 5 xV^G(~M 3 R-^G AB d A <5>d B <5> - V{$)) , 

V($) = ^($ 2 -V) 2 + A , (1) 

where X A (A = 0, 1, 2, 3, 4) is the 5D coordinates and we also use the notation 
(X A ) = (x M , y), fi = 0,1, 2, 3. X 4 = y is the extra axis which is taken to be 
a space coordinate. The signature of the 5D metric Gab is ( — h + + +)• $ 
is a 5D Higgs (scalar) field, G = detCrAB; R is the 5D Riemannian scalar 
curvature. M and V(<&) are the 5D Planck mass and the Higgs potential 
respectively. The three parameters X,vq and A in V(3>) are called here vac- 
uum parameters. A(> 0) is a coupling, t>o(> 0) is the Higgs field vacuum 
expectation value, and A is the 5D cosmological constant. 

Assuming the Poincare invariance in the brane, the line element can be 
written as 

ds 2 = e-^^dx^dx" + dy 2 = G AB dX A dX B , (2) 

where r]^ = diag(— 1, 1, 1, 1). (Gab) is explicitly written as 
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Vnv, 



(G AB ) = f Q '^' !M , v^G = e" 4CT . (3) 



The 5D Einstein equation gives us 

-6MV) 2 = --($') 2 + ^ > (4) 

3MV' = ($') 2 , (5) 

where ,u =d/dy" and we consider the case that $ depends only on the extra 
coordinate y; $ = $(?/). The above equations are 1) translation invariant 
(y — > y+ const.); 2) Z 2 symmetric (y — > —y,u' — > —0"',$ — > ±$); 3)even 
with respect to the ^-reflection ($ <-> — $). Besides they are 4) global scale 
invariant, when some vacuum parameters change appropriately; 



y-*ky,\^-rp,A->— ,v ->v , (6) 

where k is a constant. This invariance says the scale of y can be adjusted by 
the scaling of A and A. Note that the scaling power is independent of their 
mass-dimensions. | 

Eq.(^) gives an important positivity relation, 

a" = -*-r($') 2 > > 

3M 3 {a'\ y=y2 -a'\ y=yi }= (<$>') 2 dy > , y x < y 2 , (7) 

where non-singularity of o" is assumed in the region y\ < y < y 2 . This 
relation will serve as a consistency check of the solutions.^] We will also focus 
on the (non) singularity of the bulk curvature: R = —8a" + 20a' . 

As the extra space (the fifth dimension), we take the real number space 
R = (— oo, +oo). This is a simplified version of the original RS-model[l[] and 
was examined in the subsequent workpl. 



2 The mass-dimensions of the vacuum parameters, A, A and vq, are (mass) , (mass) 5 , 
and (mass) 3 / 2 respectively. 

3 For a general argument about the consistency of the domain world configuration, see 

El. 



3 Exact Solution for the No Potential Case 

Let us consider the case of no potential, V A = 0;A = A = 0. Then eq.f^J^) 
reduce to 

-6MV) 2 = --($') 2 , 3MV = ($') 2 • (8) 

cr(y) is solved as 

1 1,1-4-%! 

a =A^Ty • "=-i h 4r • B> ° • < 9 » 

where A{— oo < A < oo) and -B(> 0) are integration constants. The con- 
stant, A, comes from the translation invariance of (D). i? comes from the 
global scale invariance of (H). The line element is given by 



rfs 2 = J- — -^-ri^dx^dx" + dy 2 . (10) 

$' and $ are solved as 

*< = ± « , * = T ^,„^M , c>0 , (U) 

A — Ay 2 C 

where C(> 0) is another integration constant. The plural signs come from 
the evenness of (|J) under the " <J>- reflection" : $ <-► — $. The Higgs field, $(w), 
does not go to a constant in the asymptotic region \y\ — > oo, which should 
be compared with other solutions obtained later. 
The 5D Riemann scalar curvature is obtained as 

& = V^W <0 ' (12) 

which is 1) negative definite, 2) singular at y = A/A and 3) vanishes for 
\y\ — > oo. The metric ( |TU| ) has no horizon, hence this curvature singularity 
is a naked one. 

The obtained solution has unwanted properties and can not be used as 
the brane world model. This model is too simple. It has, however, some 
common or comparative features to the more realistic solutions of later sec- 
tions in some points such as 1) the appearance of some integration constants 



in relation to some symmetries of the field equations, 2) the plural signs, 3) 
(naked) curvature singularity, 4) no horizon. Furthermore a suggestive rela- 
tion to the R-S metric can be found by considering the region : \A\ 3> \y\. In 
this case the metric (|10D reduces to 



ds 2 « J^-e-^ri^dx^dx" + dy 2 , \A\ > \y\ , (13) 

which looks like the RS-type metric although we cannot take \y\ — > oo. 



4 Exact Solution for the Cosmological Term 
Dominated Case 

Let us examine a little more general case, that is, the cosmological term dom- 
inated case: A = 0. 

-6MV) 2 = -~($') 2 + A , (14) 

3MV = ($')' • (15) 



We have interest in the case: $ — > const, as y — ► ±oo. Then, from Q14T) , 
A <0. 

From above equations, we get 3M 3 cr"/(12M 3 cr / + 2A) = 1, which can be 
integrated as 



-=^T37^ > "Vo^ ' (16) 

where A is an integration constant (— oo < A < oo) which comes from the 
translation invariance of ( |il] , |T5| ) . 

(i) A < 

The solution (|T6|) can be written as 

a' = -cjtanh 4w(y-y*) , \A\ = e~ 8ujy * . (17) 

It is attractive that this solution is non-singular. But it contradicts with the 
positivity relation (^). Hence we conclude, in this case (i), there does not 
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exist a consistent solution. 

(ii) A = 

This case is solved as 

a' = oj , a = ujy + B , $ = v , (18) 

where B and Vq are another integration constants. The 5D curvature is a 
positive constant everywhere. 

3 M 3 v ; 

The geometry is Anti de Sitter space. 

(hi) A > 

(|i~6f ) is singular at y = — (ln.A)/8u> = y*. The solutions are obtained as 

a' = —uj coth 4u)(y — y*) , 



sinh(4c^(y - y*)) 
$ = ±^— - ln|tanh(2o;(y-y*))|H-^o , (20) 

where t>o is an integration constant. Fields asymptotically behave as a' — > 
=pjj, $ — > v o when y — y* — ► ±oo. The 5D scalar curvature is obtained as 

M 3l 3 {sinh(4cu(y-y,))} 2/ ' l ; 

which is singular at y = y*. It asymptotically behaves as i? — ► — y^ when 
|y ~~ V*\ "~ > °°. Integrating a' in eq.fl2Tf), the line element is obtained as 

ds 2 = Byf\ sinh(4w(y - y*))\r]^dx fl dx u + dy 2 , (22) 

where B is another integration constant. From this result, we see there is 
no horizon. The curvature singularity at y = y* is the naked one. In the 
region far from the singularity (Auj\y — y*| ^> 1), the line element can be 



approximated as ds 2 ~ (B / \/2)e 2 ^ y ^ y *^r] flu dx tJ 'dx l/ + dy 2 , which is similar to 
the RS metric except the sign of the exponent. 

We can check the solution of (iii) goes to (ii) as A — > +0(y* — > oo). We can 
further check the solution (iii) is continuously connected to the no potential 
case of Sec. 3 by taking A — > —0. The parameter A in Sec. 3 (= A no ~P ) is 
related to A of (0) as l/A n(> P ot = u(l + A)/(l-A). In the point that $,ct' 
and R become (non-zero) constants in the asymptotic region, the solution of 
this section approaches, compared with Sec. 3, a realistic one. 

In Sec. 3 and 4, we have obtained the exact solutions. They become RS- 
type solutions for special regions of y. ( It suggests that further "deforma- 
tion" of the potential makes us find the domain wall solution. Indeed it will 
do so. ) The solutions, however, still have a bad property of (naked) curva- 
ture singularity. To seek a non-singular solution, we must take into account 
all vacuum parameters X,v and A. That is the following subject. 

5 Domain Wall Solution 

Let us solve the 5D Einstein equations ([|,|5|) for the general case of vacuum 
parameters. We impose the following asymptotic behaviour (boundary con- 
dition) for the Higgs field $(y). 

<5>{y)^±v , y^±oo (23) 

This means $' — > 0, and from (|^), o" — > 0. From this result and (0), we are 
led to a' — ► ±u,a — > u\y\ as y — ► ±oo, where u>(> 0) is some constant. It 
can be fixed, by considering y — > ±oo in (HI), as 




uj = J M~2 , A<0 , (24) 



where we see the sign of A must be negative, that is, the 5D geometry must 
be anti de Sitter in the asymptotic regions. We may set M = 1 without am- 
biguity. (Only when it is necessary, we explicitly write down M-dependence.) 
We notice, in the results of previous sections, a' and $ behave in a com- 
parative way. Let us take the following form for a'(y) and $(|/) as a solution. 



n=0 \ AU ' -V- 



OO J 

^y) = voT,7^r I ru^H k y + l )} 2n+1 , (25) 

where c's and d's are coefficient-constants to be determined. 4 The odd- 
power terms are only taken here using the Z 2 symmetry of (f|,|5D and the 
boundary conditions explained above. The free parameter I comes from the 
translation invariance of (f|) and (|5|). A new mass scale k(> 0) is introduced 
here as a free parameter to make the quantity ky dimensionless. The freedom 
comes from the global scale invariance of @ and (^). The physical meaning 
of 1/k is the thickness of the domain wall. The parameter k plays a central 
role in the dimensional reduction scenario 0. The distortion of 5D space- 
time geometry by the existence of the domain wall should be small so that 
the quantum effect of 5D gravity can be ignored and the present classical 
analysis is valid. This requires the condition [[jj 

k < M . (26) 

Besides M — 1, we can also take k = 1 without ambiguity (keeping the 
relation (pS| ) in mind). The coefficient-constants c's and d's have the following 
constraints, 



oc 



C2n+1 -, v^ "2n+l 



"n^ + l)! ' "^(2n + l)! ' l ' } 

which are obtained by considering the asymptotic behaviours y — > ±oo in 

We first obtain the recursion relations between the expansion coefficients, 
from the field equations (|4|) and (||). For n > 2, they are given by [|7j 

C2n+1 C2n - X V ° 2 -{D' n -2D' n _ l+ D' n _ 2 ) , 



(2n)! (2n-2)! /^X/2 



Cn-l : 




•^-(D / „-2 J D / „_i + J DV 


-2) + -T-r-(E n -2 - 2D n _i) , 


where 








n 
m=0 


(2r 


«2n-2m+l«2m+l 

i-2m+l)!(2m + l)! ' 


p./ y^ "2n-2m+l«2m+l 

n "^ (2n-2m)!(2m)! 



(28) 



4 Normalization of C2 n +i is different from that of ref.H. The relation is C2 n +i — 
(k/uj) x c 2 „+i of ref.0. 



" r r " 

rt _ V^ ^2n-2m+1^2m+l p _ \^ p, p, /oq\ 

^ (2n-2m + l)!(2m+l)! ^ 

The first few terms, (ci, di), (c 3 , rf 3 ), are explicitly given as 

Xv,\ 




d 1 = ± — \/A+— -- , d = — ^^(A + 
' "' 3,/-A/6 



d f = 2 + {|(A + ^f) - W} , ^ = 2 + {^(A + *f)- 2AV} , (30) 
d\ 3 4 ci 3 4 

where ± sign in d\ reflects $ <-» — $ symmetry in (B) and (|5]). We take 
the positive one in the following. We can confirm that the above relations, 
( |28| ) and (|30"D, determine all c's and gPs recursively in the order of increasing 
n. In (|28|) and (|30"D , M = fc = 1 is taken for simplicity. Their dependence 
is easily recovered by A — > A/k 2 M 3 , A — > XM 3 /k 2 , v — > fo/v 7 ^". Note 
that all coefficients, derived above, are solved and are described by the three 
(dimensionless) vacuum parameters. Among the 3 parameters, there exist 



2 constraints from (|27j). Hence the present solution is one-parameter family 
solution. 

In order for this solution to make sense, as seen from the expression for 
d\, the 5D cosmological term A should be bounded both from below and from 
above. 

-^<A<0 . (31) 

The presence of this relation says the non-singular solution presented here 
cannot continuously connect with the singular solutions of Sec. 3 (A = A = 0) 
and of Sec.4(A = 0). 



6 Evaluation of Coefficients and Numerical 
Check of Analytic Results 

We present here the results of concrete values of c's and d's for two input 
values vq = 1.0 and v = 1.6. We solve constraints (P7|) by taking the 
first 7 terms (up to n=6th order). The most important point is to confirm 
the convergence of the infinite series Z^o (2n+i)p anc ^ S^lo (2n+i)\ > wm ch 
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Fig.l The values of c 2n+ i/(2n + l)!(blob),d 2n+ i/(2n + 1)! (small circle). The 
large circles show (0.25) n . ( n — 0, 1, • • • , 6 .) v = l.O(input). 



guarantees the present boundary condition. The 6th-order approximation 
calculation gives the vacuum parameters as 



v = 1.0 (input) 
Vq = 1.6(input) , A 



A = -0.272 , A = 2.88 
-1.50825 , A = 1.49925 



(32) 



[] The input value v = 1.6 is quite near the case of [0], and the above results 
of A and A are consistent with the previous results. Both cases give similar 
behaviours, hence we present further results only for vq = 1.0. The obtained 
values of the coefficients are shown in Fig.l. In the figure we also plot the data 
from the geometrical series: 1/(1 — x) — 1 +x + ■ • - at x — 0.25. Comparing 
them, we can recognize the (rapid) convergence of the coefficient-series (up to 
this approximation order). Note that two series, { (^T+iv. } anc ^ { (2n+iy. }' are 
'oscillating'. Using these results, the analytical results of the scalar field $(y) 
and the warp factor <j'(y), (|25|) , are shown in Fig.2. f\ The (5D) Riemann 






C2n + 1 



(i-ELo- 



< 1.1 x 



The constraint is satisfied as (1 - ^ n=0 {2n+iy l t^- z^„=o (2 n+i)P 
10 -7 (^o = 1.0), 1.6 x 10 -9 (vq = 1.6). In the search of the solution (A, A), the mesh-size 
of the parameter-space determines the "precision". The high precision for the vo = 1.6 
case shows the mesh-size is taken much smaller than that of vq = 1-0. 



We use the truncated version of the expressions (E5|) by the first seven terms 
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k y 



Fig. 2 The analytic result of a' /u> (dashed line) and $/t>o (normal line). Both 

are odd with respect to y <-► —y. The graphs are depicted by using ( f23| ) 
with the 6-th order approximation, vq = l.O(input). The horizontal axis is 

ky. 



scalar curvature is also shown in Fig. 3. It shows the curvature is negative 
inside the wall, whereas positive outside. There exist two points, ky ~ ±1, 
where the curvature vanishes. We see the present solution is non-singular 
everywhere. The presence of a dip around y = (Fig. 3) clearly says that the 
domain wall exists there. 

The kink b.c. for the Higgs field ( p3|) guarantees the stability of the 
present solution. It reflects the b.c. for a' as specified by the parameter 



value (|24]). We can also see the stability from the behaviour of R (Fig. 3) as 
follows. From the expression (|25| ) and the b.c. (|27D , R has the following b.c. 
in the IR region. 



\a \ —>■ uj , a —* , — - ^ 20 



as \y\ — » oo 



We can also see, from (pp]), the b.c. in the UV region is 



(33) 



a 



a 



ujkci 



R 



ijj- 



8k 
Cl 



as 



2/^0 



(34) 
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k y 



Fig. 3 (5D) Riemann scalar curvature R/u 2 in the 6th order approximation. 
It even with respect to y <-> —y. The horizontal axis is ky. Vq = l.O(input). 



where c\ > from the b.c. of a'. Due to the continuity, R must have a dip 
around the origin with a finite thickness. [] 

As the coupled differential equations for $(?/) and cr'(y), the equations 
(^|||), have the standard form of the numerical analysis, that is, Runge- 
Kutta method. We can numerically solve them without any ansatz about the 
form of the solution. In this numerical analysis, the following two points are 
important: 1) the choice of three parameters v o, A and A; 2) the choice of the 
initial conditions, o~'(y = 0) and $(y = 0). As for the point 1) we can borrow 
the values obtained in the 6th order approximation (|32|). As for 2), due to 
the required Z 2 symmetry (the oddness under y — > —y) for $(?/) and o~'(y) 
(|25|), we can take y = as the initial point of y and the initial conditions 
$(0) = o~'(0) = 0.0. The numerical result is shown in Fig. 4. It shows the 
analytic solution of Fig.2, based on Sec. 5, is reproduced very well. The 
numerical output data stop at ky ~ 3.0 with producing imaginary values. 
This occurs because keeping the positivity, $'(?/) > 0, in the numerical 
analysis becomes so stringent in the infrared region. The quantity becomes 
so small in that region and vanishes at k\y\ = oo. We understand that further 

7 For the input v = 1.0, the obtained values, A/k 2 M 3 = -0.272, ci = 1.4, give 
—8kci/uj = —52.6 which is shown in Fig. 3. 
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ky 



Fig.4 The numerical results for cr'/a;(up) and $/t>o(down). They are 

obtained by Runge-Kutta method. One step value along %-axis is 0.05. 

About 65 points are plotted for each line in the figure. The initial point is 

y = 0. The horizontal axis is ky. 
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higher order calculation is required for the values of A and A (for an input 
value vq) in order to extend the valid region furthermore. 

7 Discussion and conclusion 

The assumption about the convergence of the series (|25|), which was assumed 
in 0] , is strongly confirmed. The present result of the 6th order calculation 
does not so much deviate from that of the 2nd order one in . It says the 
truncation approximation of (|27| ) is valid. 

We point out some results which are potentially important in phenomenol- 



ogy. The cosmological term has both the upper and the lower bound (31). 
It is expected to be useful when we fix the parameters A, t>o and A. The sign 
change of the curvature near the "surface" k\y\ ~ 1 (Fig. 3) could become an 
important check point of the confirmation of this model. 



From the point of singularity resolution^], the curvature singularity ap- 
pearing in Sec. 3 and Sec. 4 is resolved in Sec. 5 by a sort of "deformation" 
(adding the potential terms appropriately). In the procedure two constraints 
appear among the three vacuum parameters. 

The present standpoint is that the domain configuration should be re- 
alized in the non-singular geometry. The approach based on the singular 
geometry, such as the original one0], can be regarded as a temporary stage 
of the development. The singularity, often expressed by the (^-function intro- 
duced by hand, is expected to be derived by some definite limiting procedure, 
say, the thickness goes to zero: k~ x — > 0. If the string or the D-brane is the 
fundamental constituents of nature, such extended objects are strongly ex- 
pected to behave smoothly in the UV-region. Because the domain world 
physics can be regarded as a transitive approach from the field theory to the 
string-brane theory, we believe seeking non-singular solutions is important in 
the development. 
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